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ABSTRACT 

We investigate the red supergiant problem, the apparent dearth of Type IIP super- 
nova progenitors with masses between 16 and 30 Mq. Although red supergiants with 
masses in this range have been observed, none have been identified as progenitors in 
pre-explosion images. We show that, by failing to take into account the additional ex- 
tinction resulting from the dust produced in the red supergiant winds, the luminosity 
of the most massive red supergiants at the end of their lives is underestimated. We 
re-estimate the initial masses of all Type IIP progenitors for which observations exist 
and analyse the resulting population. We find that the most likely maximum mass for 
a Type IIP progenitor is 21^1 Mq. This is in closer agreement with the limit predicted 
from single star evolution models. 
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1 INTRODUCTION 

Nothing in the life of a massive star becomes it like the 
leaving it. A supernova (SN) is one of the most impressive 
spectacles that the Universe can afford an astronomer. How- 
ever, there is some uncertainty as to the range of stars that 
undergo this most spectacular demise. Observations indi- 
cate that stars evolve into red supergiants if their initial 
masses are between about 12 and 30 Mq l|Levesaue et al.l 
I2OO6I ). Stars that are less massive undergo second dredge-up 
and end their lives as asymptotic giant branch (AGB) stars 
l|Eldridge et al.ll2007l ). Mass-loss rates increase with mass 
and stars that are more massive suffer enough mass-loss to 
remove their hydrogen envelopes before they die. They be- 
come Wolf-Raye t stars, naked helium stars with thick winds 
l|Crowtheij|2007h . These explode as hydrogen-free Type Ib/c 
SNe, rather than the more common Type II SNe that are be- 
lieved to result from the d eath of red supergiants (|Filippenkd 
ll997l : ISmartt et aljlioosh . Type II SNe are in turn divided 
into Types IIP, IIL, Iln and lib. Types IIP and IIL are iden- 
tified by their SN light curves. The former have a plateau 
and the latter show only a linear decline. The plateau is 
the result of the photosphere maintaining a cons tant radius, 
movi ng inward in mass as the ejecta expands (|Filippenkd 
Il997l h This in turn is due to the ionised hydrogen recombin- 
ing. Type lib SNe have weak hydrogen lines and light curves 
similar to hydrogen-free Type lb SNe, implying that they 
contain only a small percentage of hydrogen. There seems 
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to be a sequence from Type IIP to Type IIL to Type lib 
SNe, driven by increased mass-loss and a consequently re- 
duced mass of hydrogen in the ejecta. Type Iln SNe are dis- 
tinguished by narrow line hydrogen emission, the result of 
shock interaction with circumstellar material (|Smith et al.l 
[ioil). The mass limits have some dependence on metal- 
licity because metal-rich stars have higher mass-loss rates 
and thus the minimum initial mass of W olf-Rayet stars is 
appro ximately 25 Mq at Solar metallicity (jEldridge fc Tou^ 
12004 ). We therefore expect Type II progenitors up to this 
limit. 

The red super giant problem was first reported by 
ISmartt et all l|2009l ). They compared 20 Type IIP SN pro- 
genitor detections and non-detections with stellar evolu- 
tion models to determine the minimum and maximum ini- 
tial mass limits for the progenitors of these SNe. They 
found that the minimum mass required for stars to ex- 
plode as Type IIP SNe was S.Slj 5 Mq, which is consistent 
with the observed upper lim i t for white dwarf for mation 
f Weidem ann fc Koeste J 1 19821 : IWiUiams et~ai1 |200i). More 
surprisingly, they found an upper limit of IGIJ^ ^Mq, a 95 
per cent confidence upper limit of 21 Mq. In essence, the red 
supergiant problem is that this estimate is well below the 
maximum mass estimated for red supergiants. There there- 
fore appears to be an absence of higher-mass red supergiant 
SN progenitors, leaving the fate of stars with masses between 
16 and 25-30 Mq uncertain. 

It should be mentioned that the alternative to deduc- 
ing masses from progenitor models is to model the su- 
pernova directly. This has the advantage of not requiring 
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pre-explosion images, although it does depend on detailed 
follow-up observations of the supernova itself. Utrobin and 
Chugai have used a one-dimensional h ydrodynamic code 
llUtrobinll2003 ') to model the SNe 2005cs jutrobin fc Chugail 
I2OO8I ) and 2004et (|Utrobin fc Chugail |2009| ). Both have de- 
tected progenitors and in this paper we deduce initial masses 
of 9^4 and 12^} Mq respectively. This compares with 18 
and 28 Mq for the hydrodynamic masses. If these masses 
are correct, the red supergiant problem ceases to be. How- 
ever, three earlier independent attempts to model the pro- 
genit or of 2005c s gave masses of 9^2 (jMaund et al.l 
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(iLi et al.|[2006l ') and between 6 and SMg 



I Eldridge et alj|2007h . a disc repancy noted by t he au thors. 
There are other approaches; iDessart fc Hillien l|201ll ) have 
constructed models that include the nebula phase of the su- 
pernovae, which is more amenable to detailed simulation. 
These SNe modelling approaches have great potential but 
wait on more sophisticated models. 

Another approach is to consider the ratios of the differ- 
ent types of SNe. One then integrates the initial mass func- 
tion (IM F) to find the limits that provide the desired number 
of stars. ISmith et al.l l|201ll ) considered the fractions of core 
collapse SNe from the Lick Observatory Supernova Search 
(LOSS) and observed that the proportion of Type Ib/c SNe 
was much too high for all of them to be the result of sin- 
gle star evolution. In addition, with the standard Salpeter 
IMF, the observed fraction of Type IIP SNe was such that it 
could be produced by single stars with initial masses in the 
range 8.5 — IS.TMq. These facts imply that binary interac- 
tion allows the production of hydrogen-free p rogenitors at 
lower masses than would otherwise be the case. ISmith et all 
suggest that if binaries are included the upper limit 
for red supergiants would be in the range of 18 to 24 Mq. 

However, this study found a lower fraction of Type IIP 
SNe than previous work. The core collapse SNe broke down 
as 48 per cent Type IIP and 22 per cent Type Ib/ c, com- 
pared with 59 and 29 respectively for the survey of 'Smart t| 
Hooi). The reason for this is not clear. iSma rtt (2009) con- 
sidered a 28 Mpc volume-limited survey using all de tected 
SNe within that volume, whereas ISmith et al. (|201ll ) used 
60 Mpc and only those SNe detected in LOSS. Type IIP SNe 
tend to be dimmer than other types so those at large dis- 
tances may have been missed. While these selection effects 
were accounted for, such adjustments are, by their nature, 
very uncertain. Another reason for the difference could be 
that ISmith et all l|201lh had more complete spectroscopic 
and light curve follow-ups and so had greater accuracy in 
their classifications. 

It is possible to obtain the required SNe fractions with 
an ap propriately chosen stellar population. lEldridge et al.l 
l|201ll ) showed that a population composed of a mix- 
ture of binary and single stars could e xplain the rates of 
ISmartt et al] (|2009l ). [Smith et al.l (|201lh used their own bi- 
nary population models to explain their results. The uncer- 
tainity in the rates means that it is hard to use these frac- 
tions to tightly constrain the progenitor mass range, partic- 
ularly when the shape of the IMF means that small changes 
in the lower limit result in large changes in the population 
fractions. In contrast large changes in the upper limit re- 
sult in small changes in the population fractions. Hopefully 
future surveys will resolve the discrepancy. 

If we accept the existence of the red supergiant prob- 



lem then we must consider a number of possible explana- 
tions. First, that these massive red sup ergiants form blac k 
holes with faint or non-existent SNe (e.g. lHeger et al.ll2003l ). 
Secondly, that their envelopes a re unstable and eject a 
large amount of mass pre-SN (e.g. lYoon fc Cantielloll201ol ). 
Thirdly, that they explode as a different type of SN (e.g. 
iKotak fc Viii^booel ). We suggest a fourth explanation, that 
the mass estimates may be systematically inaccurate at 
the high-mass end. Mass estimates are based on mass- 
luminosity relations, so extra intrinsic extinction close to the 
red supergiant progenitors would give reduced luminosities 
and hence lower predicted masses. ISmartt et al] t2009,) were 
careful to provide extinction estimates when possible from 
measurements of nearby supergiants and from the supernova 
itself. These could be underestimates. If red supergiants pro- 
duce extra dust local to the star, it would be destroyed in 
the supernova explosion. This is very plausibl e. It is known 
that red supergiants form dust in their winds (|Massev et al.l 
l2005h . Furthermore, IR interferometry has sh own that this 
dust c an be found very cl ose to the star itself (|Danchi et al.l 
ll994h . lSmith et all l|2011^ also suggest dust as a solution to 
the red supergiant problem. 

In this work, we first describe the theoretical models 
from which we derived mass-magnitude relations. We then 
show how these were used to deduce masses from a popu- 
lation of Type IIP SN progenitors. Finally, we deduce the 
most probable upper and lower mass limits and present our 
conclusions. 



2 SIMULATING RED SUPERGIANTS 

We begin by considering the Cambridge STARS code, the 
source of our stellar models. These models were processed 
to generate obs erved colours with t he BaSeL V3.1 model at- 
mosphere grid jWestera et al.|[2002h and the relevant broad- 
band filter functions. We used the dust production rate ob- 
served by iMassev et al.l |200^ to estimate the amount of 
circumstellar extinction that would then manifest. This al- 
lowed the calculation of mass-magnitude relations both with 
and without the inclusion of circumstellar dust. We also con- 
sider the nature of dust production, including non-spherical 
behaviour. 



2.1 The Cambridge STARS code 

The Cambridge STARS co de was originall y developed by 
Peter Eggleton in the 1960s (lEggletonlligTll ). It uses a non- 
Lagrangian mesh, where the mesh function ensures that the 
points are distributed so that no quantity of physical inter- 
est is allowed to vary by a large amount in the intervals. 
The code has been gradually improved and updated and 
the work in this paper is b ased on the code described by 
IStancliffe fc Eldridgd l|2009l ) and those referenced by them. 

Con vection is included in the code by the mixing length 
theory of lBohm-Vitensj (|l958l ). with a solar-calibrated mix- 
ing length parameter of a = 2.0. Convective overs hooting is 
obtained with the method of ISchrder et al.l (|l997l ). with an 
overshooting parameter of Sqv = 0.12. This method involves 
the addition of a 5 term to the adiabatic gradient, allowing 
mixing to occur in regions that are weakly stable by the 
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Schwarzchild criterion. The code follows the chemical evolu- 
tion of ^H, ^He, ^He, ^^C, "N, ^^O and ^''Ne, together with 
structural variables. 

We use the m ass-loss scheme described by 
lEldridge fc ToutI (|2004l '). For main-sequence O B stars the 
mass- loss rates are calculated according to IVink et al.l 
and for all other stars we use the rates of 
Ide Jager et al.l (|l988l '). where the metallicity scaling goes as 
(z/zq)"'^. This theory is older but the the rates have been 
recently tested for red supergiants and been shown to be 
still the best rates available for them (jMauron fc JosselinI 

We have created a library of evolution models, with val- 
ues of Z, the metallicity fraction by mass, equal to 0.02, 0.01, 
0.008 and 0.006. The se cover the m e tallici ty range of the IIP 
progenitors given bv ISmartt eTsIl (|2009l ). The fractions of 
hydrogen and helium were determined on the assumption of 
constant helium enrichment from the primordial condition 
of X=0.75, y=0.25 and calibrating to a Solar composition 
of J£:=0.70, y=0.28 and Z=0.02, i.e. that Y = 0.25 + 1.5Z. 
We evolve our models through to core neon burning, a few 
years before core collapse. 



2.2 The colour— magnitude diagrams 

The Cambridge STARS code outputs the physical parame- 
ters for a stellar model at each time-step. These include the 
bolometric luminosity and the surface temperature. Because 
the stellar observations are in the form of colours and mag- 
nitudes, either the luminosity and the surface temperature 
must be estimated from the observations or a method of cal- 
culating colours and magnitudes from the models must be 
developed. We choo se the latter c o urse b ecause this requires 
fewer assumptions. ISmartt et al.l (|2009l ') used the opposite 
approach and, for that reason, the masses they deduced 
sometimes differ by a small amount. W e use the method de- 
scribed by lEldridge fc Sta nwavl (|2009l ) to process the code 
output and calculate magnitudes to compare to those ob- 
served. The BaSel v3.1 grid of model atmospheres is ar- 
ranged over surface temperature and effective gravity. Us- 
ing the values of these variables from the code, we obtained 
appropriate templates for the SEDs for each model at each 
time-step. We then applied the filter functions to extract the 
magnitudes in the various bands. This allowed us to produce 
colour-magnitude diagrams from the evolution tracks. 

To check the validity of our synthetic colours we com- 
pared our models to red supergiants ob served in the Mag- 
ellanic Clouds by Levesgue et all (|2006l ^ and in the Milky 
Way bv I Levesgue et all ( 2005f ). These observations included 
estimates for the total extinction based on spectrophoto- 
metric modelling, so it was possible to process the data to 
get the absolute colours and magnitudes in the absence of 
extinction. This meant that the models did not yet have 
to take into account the effects of circumstellar extinction. 
This is shown in Figure [T] We see that the models agree 
well with the observations, with the red supergiants appear- 
ing towards the end of the evolution tracks and ranging in 
mass between about 10 and 30 Mq . Some of the SMC stars 
are cooler than predicted and may have higher metallici- 
ties. This is not unexpected with a large and heterogeneous 
population. 

To test the validity of the models further, we took the 





Figure 1. Evolution tracks for every integer mass between 5 and 
30 Mq , with the multiples of 5 indicated by dashed lines. The 
crosses are the observed red supergiants. The SMC stars are at 
the top and the models use Z = 0.004. The LMC stars are in the 
middle and the models use Z = 0.01. The Milky Way stars are at 
the bottom and the models use Z = 0.02. 



data and made cumulative frequency plots in Mk ■ We made 
similar plots from the models by first identifying which stars 
became red supergiants during their lifetimes. We required 
that V — K > 2.5 so as to get the reddest stars. The lower 
limit in K was set to be —9.5 for the Clouds models and 
—8.5 for the Galactic models. These limits are reasonable for 
red supergiants and were chosen to reflect the distribution 
of the observations. We weighted the K magnitudes by the 
timestep to reflect a greater probability of observation, and 
the Salpeter IMF. We used metallicities of Z = 0.004 for the 



© 0000 RAS, MNRAS 000, [l]-?? 



4 Joseph J. Walmswell, John J. Eldridge 



SMC, Z = 0.01 for the LMC and Z = 0.02 for the Galactic 
population. The comparison is shown in Figure [2] 

The agreement is quite good, although there are fewer 
very luminous stars in the LMC and SMC sets than this sim- 
ple application of our models predicts. This is probably be- 
cause all three data sets were not compiled to reflect a stellar 
population but as an observationally convenient group of red 
supergiants. In addition, because we weight by the timestep 
and the IMF only, we are implicitly assuming a constant 
rate of star formation, a doubtful supposition. We have only 
considered single stars and binary int eractions may change 
these predicted synthetic frequencies (|Eldridge et al.l [20081 ) 
and we may be underestimatin g the mass-loss ra.tes of the 
most luminous red supergiants (jYoon &: Cantielldlioioh . 

2.3 Circumstellar dust 

When the surface temperature of a red supergiant falls below 
about 5000 K, dust begins to condense out of t he stellar wind 
at a di stance of around 5 — lORstar ~ 1000-Rq (|Massev et al.l 
l2009al ). It might be expected that the amount of dust pro- 
duction would correlate with the mass-loss, which in turn 
correlates with the luminosity , because thi s is res ponsible 
for th e stellar wind (jvan Loon et al.l |2005| ). iMassev et al.l 
l|2005ll showed that the dust production rate correlates 
with the bolometric luminosity, with a least squares fit of 
logio(Mdust) = — 0.43Mboi — 12.0, where the dust produc- 
tion rate has units of M0/year. 

The dynamics of the wind and the dust is complicated 
and simulat i ons i ndicate an absence of spherical symme- 
try IWoitkd i 20061 ) found that various instabilities such as 
Rayleigh-Taylor or Kelvin-Helmholtz lead to the formation 
of arcs and caps of dust despite the spherical initial condi- 
tions. This means that one would expect variation in the 
observed mass-loss and extinction that is entirely due to 
the behaviour of the dusty wind along the line of sight. This 
cannot be accounted for with the STARS code, because in 
the absence of observations of the dust, one can only use the 
relation between the luminosity and the dust production to 
estimate an average extinction. The additional variation due 
to the lack of spherical symmetry means an additional source 
of uncertainty. 

The amount of extinction that is due to circumstellar 
dust depends on the past history of the star, which has a 
continuous loss of mass over time. In contrast, the STARS 
code calculates the properties of the stars at intervals deter- 
mined by the time-step. To account for this, we calculated 
the launch velocity of the dust and the dust-production rate 
for each time in the output data. This, together with the 
stellar radii, when interpolated, gave the distri bution of dust 
with distance from the st ar. Following Massev et al.l (|2005| ) 
we begin by referring to IWhittetl (|2003l )" who showed how 
the extinction owing a thin shell of dust can be determined. 
If one assumes a dust grain density of s = 2500 kg m~^, 
applicable to low-density silicates and a refractive index of 
m — 1.50, one can obtain the extinction Ay in terms of the 
path length L. 

«M3,7xlO-„=!±|f. (1, 

After the substitutions are made, one obtains the den- 
sity in terms of the path length and the extinction only. 




Figure 2. The cumulative frequency diagrams, in K band mag- 
nitude, for red supergiants observed in tlie SMC, the LMC and 
the Milky way. They are compared with cumulative frequency 
curves from synthetic populations derived from the same models 
used for Figure [T] 



Pd = (3.1 X 10-*)^. (2) 

If a thin shell of dust is used, the thickness of the shell 
5R cancels with the path length. The extinction caused by 
the shell is then in terms of the dust mass Md and the radius 
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Table 1. Extinction at standard and HST pa ss-bands relative to 
tlie V band, as given bv lCardelli et alj lll989h . 



Filter 


Aa/Av 


U 


1.569 


B 


1.337 


V 


1.000 


R 


0.751 


I 


0.479 


J 


0.282 


H 


0.190 


K 


0.114 


F555W 


0.996 


F606W 


0.885 


F814W 


0.597 



Av = 



(3.2 X 10^)LMd 
(3.2 X 10=*)Md 



(3) 

The dust was modelled as a series of these thin shells, 
over which the total extinction was integrated at each 
timestep. The Ay was then used to calculate the extinction 
in the other pass-bands by using the extin ction law and as- 
sociated ratios described in lCardelli et al.| (1989). These are 
shown in Table [1] 

The density of dust falls off according to an inverse 
square law as each shell moves outwards in the stellar wind. 
Therefore only material very near to the star has a signif- 
icant effect. We find about 95 per cent of the extinction is 
due to material closer than 50Rstar in all the models. This 
means that the dust is unlikely to affect the SN because it 
is rapidly swept up and destroyed in the explosion. The in- 
verse square law also means that the distance at which dust 
first forms in the wind is important. We chose to set this to 
10 Ratar, au Upper estimate, to ensure that our extinctions 
are modest underestimates. 



3 THE SUPERNOVA PROGENITORS 

W e use the comp i lation of SN detections and non-detections 
of ISmartt etaU Hooi). All progenitor information can be 
found in that paper and the re ferences therein. We supple- 
ment this with SN 2009md (Er aser et al.ll2011al '). 

The metallicities are based on neighbouring 0/H 
number ratios, where [0/H] = logio(0/H) -I- 12. For 
[O/H] > 8.7, the Z = 0.02 models were used. Similarly, 
for 8.5 < [O/H] 8.7, Z = 0.01, for 8.4 < [O/H] s: 8.5, 
Z = 0.008 and for 8.2 < [O/H] sC 8.4, Z = 0.006. This is 
a more precise divi sion by metallicity than was possible for 
ISmartt et all (|2009l ). Errors are given, when known, and the 
errors for the absolute magnitude were obtained by combin- 
ing the other errors in quadrature. 

We assigned masses to the progenitors by comparing the 
absolute magnitudes, together with their error bars, with a 
plot of red supergiant final luminosities against initial mass. 
We decided to plot the minimum and maximum magnitudes 
in the lifetime of our theoretical red supergiants, the lifetime 
being between the end of core helium burning and the ter- 
mination of the model. This gave a range in luminosity over 



which a star of a given initial mass might explode. Masses 
were first obtained from the original models and again after 
processing to include the effects of circumstellar dust. The 
error in the magnitudes gave the error in the masses. 

Figure [3] shows how the predicted magnitudes for the 
Z = 0.02 models varies with initial mass in the V, R, I, and 
K bands. As expected, dust is much more of a problem at 
shorter wavelengths, when extinction is more severe, and at 
higher masses, when mass-loss is greater. 

However, although extinction is less important at longer 
wavelengths, changes in stellar luminosity have more of an 
effect. When a red supergiant becomes more luminous, most 
of this increase in output is at longer wavelengths. One can 
see from Figure [3] that the difference between the minimum 
and maximum magnitudes for a given mass is highest in 
the K band. This is particularly acute for the lower masses 
because they end core helium burning as blue supergiants, 
setting a very low minimum magnitude in the infrared. 

The lower-mass models also have higher maximum 
magnitudes than the extrapolation of the behaviour of the 
more massive stars might predict. These stars undergo sec- 
ond dredge-up, bec oming more luminou s asymptotic giant 
branch (AGB) stars l|Eldridge et al.ll200'7l 'l. In both cases, in- 
creasing the magnitude range also increases the uncertainty 
in the progenitor mass-luminosity relation. Finally the sharp 
increase in the maximum V magnitude at the high mass end 
indicates incipient Wolf-Rayet star formation. 



3.1 The non— detections 

For the non-detections, no progenitor was confidently identi- 
fied in the pre-images. This still provides useful information 
because it sets a limit on the magnitude of the progenitor. 
If the star were brighter, it would have been detected. If the 
errors in the magnitudes are approximately normal, there is 
an 84 per cent probability that the magnitude of the non- 
detected progenitor is less than the upper error bar. This is 
a sufficiently high confidence level that we took the upper 
mass limit to be the lowest mass with a magnitude range 
entirely brighter than this upper error bar. 

Some of the SNe had non-detections in several pass- 
bands. The magnitude limits quoted in Table [5] are those 
which gave the lowest upper mass limit. The other pass- 
bands merely set a higher upper mass limit which added 
no additional information. If we compare the mass limits 
deduced from the dustless models with those with dust we 
can see that the difference between the two increases with 
mass. Most notably that for SN 2003ie changes from 22 to 
25 /Mq. 



3.2 The detected progenitors 

The actual detections are fewer in number than the non- 
detections. For progenitors observed in multiple bands, a 
mass for each band was calculated and these were aver- 
aged. For SN 2009bk [Mattila et al.l l|2008ll used the progeni- 
tor SED to deduce a total extinction of Ay = 1.0 ± 0.5 and 
the absolute magnitudes take this into account. Because this 
includes any extinction from circumstellar dust, there was 
no need to use the dusty models and the predicted mass is 
the same in both cases. 
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Figure 3. The final magnitudes in V , R, I, and K from our Z=0.02 stellar models. Red indicates the dust free models and black the 
models with dust. The thinner lines are for the minimum magnitudes and the thicker are for the maximum magnitudes. 

Table 2. The observed parameters and estimated upper mass limits for Type IIP supernova progenitors that were not detected in pre- 
explosion images. We include the masses without considering extinction due to intrinsic dust, M^ustlessi E'nd the masses taking intrinsic 
dust into account M^ust- Note that we consider 2004A to be a non— detection. The observation of a progenitor was doubtful. 



SN 


Metallicity 


Distance 


Apparent 


Absolute 


Mdustlcss 


Mdust 




/dex 


/Mpc 


magnitude 


magnitude 


/Mq 


/M© 


1999an 


8.3 


18.5±1.5 


mp606W > 24.7 ± 0.2 


Mf606W > -7 ±0.3 


18 


21 


1999br 


8.4 


14.1±2.6 


mp606W > 24.91 


Mf606W > -5.89 ± 0.4 


11 


12 


1999em 


8.6 


11.7±1.0 


mi > 22.0 


Mi > -8.5 ±0.2 


18 


19 


1999gi 


8.6 


10.0±0.8 


mpeoew > 24.9 ± 0.2 


Mf606W > -5.7 ± 0.3 


12 


13 


2001du 


8.5 


18.2±1.4 


niP8i4w > 24.25 


Mf8i4w > -7.4 ± 0.2 


13 


13 


2003ie 


8.5 


15.5±1.2 


mn > 22.65 


Mr > -8.3 ±0.2 


22 


25 


2004A 


8.3 


20.3±3.4 


mF8i4W > 24.25 


Mf8i4w > -7.4 ±0.2 


14 


15 


2004dg 


8.5 


20.0±2.6 


mF8l4W > 25.0 


Mf8i4w > -6.9 ±0.3 


11 


11 


2006bc 


8.5 


14.7±2.6 


niF8l4W > 24.45 


Mf8i4W > -6.8 ±0.5 


11 


12 


2006my 


8.7 


22.3±2.6 


niF8i4W > 24.8 


Mf814W > -7.0 ±0.2 


11 


11 


2OO60V 


8.9 


12.6±2.4 


mF8i4W > 24.2 


Mf814W > -6.3 ±0.4 


10 


11 


2007aa 


8.4 


20.5±2.6 


mF8l4W > 24.44 


Mf814W > -7.2 ±0.3 


12 


13 



3.3 Other deduced masses 



3.4 The maximum likelihood limits 



Table g] lists the properties of SNe 2004am and 2004dj, 
both of which are Type s IIP but lack detected progeni- 
tors. ISmartt et al.l l|2009l ) describe how population synthe- 
sis codes were used to deduce progenitor masses from their 
parent clusters. These SNe were part of the survey and are 
included for completeness. 



The masses for the detections and non-detections have been 
drawn from a distribution of Type IIP progenitors with un- 
known parameters. If we assume that the progenitors are 
drawn from a population described by the initial mass func- 
tion (IMF), the nature of the IMF and the range of masses 
which explode as Type IIP S Ne are the importan t parame- 
ters. Following the method of lSmartt et al.l l|2009t ). we used 
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Table 3. The observed parameters and estimated masses for Type IIP supernovae that were detected in pre-explosion images. We 
include the masses without considering extinction by intrinsic dust, Mjusticssi ^'nd the masses taking intrinsic dust into account, M^ust- 



Supernova 


Metallicity 


Distance 


Apparent 


Absolute 


Mdustloss 


Mduat 




/dex 


/Mpc 


magnitude 


magnitude 


/M© 


/Mq 


1999ev 


8.5 


15.14±2.6 


mpsssw = 24.64 ± 0.17 


Mp555w = -6.7 ± 0.4 


18^3 


20+^ 


2003gd 


8.4 


9.3±1.8 


mv =25.8 ±0.15 
mi = 23.13 ±0.13 


Mv = -4.47 ±0.5 
Ml = -6.92 ±0.4 


8l? 


8tl 


2004et 


8.3 


5.9±0.4 


mi = 22.06 ±0.12 


Mi = -7.4 ±0.2 


11+1 


12+1 


2005CS 


8.7 


8.4±1.0 


mi = 23.48 ±0.22 


Mi = -6.3 ±0.3 


9+1 


9+1 

"-4 


2008bk 


8.4 


3.9±0.5 


mi = 22.20 ±0.19 
mn = 18.78 ±0.11 
mx = 18.34 ±0.07 
mj = 19.50 ±0.06 


Mi = -7.2 ±0.4 
Mh = -9.4 ±0.3 
Mk = -9.7 ±0.3 
Mj = -8.7±0.3 




121:1 


2009md 


9.0 


21.3±2.2 


my = 27.32 ± 0.15 
mi = 24.89 ± 0.08 


Mv = -4.630-3 4 
Mi = -6.92l°-^ 







Table 4. The observed parameters and estimated masses for Type IIP supernovae that have mass estimates derived from observations 
of their host clusters. 



SN 


Metallicity 


Distance 


Apparent 


Absolute 


Mnodust 


Mduat 




/dex 


/Mpc 


magnitude 


magnitude 


/Mo 


/Mq 


2004am 


8.7 


3.7±0.3 


n/a 


n/a 


12+^ 


n/a 


2004dj 


8.4 


3.3±0.3 


n/a 


n/a 




n/a 



maximum likelihood theory to find parameters that gave the 
greatest probability of generating the data. If Pi is the prob- 
ability of the ith detection or non-detection being made, 
the likelihood jC is the probability of observing the whole 
dataset. 



(5) 



Taking the natural logarithm converts the product into 
a sum and simplifies matters because maximising log^, £ is 
equivalent to maximising £. 



i=N 

^loge[Pi(m) 



(6) 



For the non-detections, the probability of the event is 
the probability that a randomly chosen star has a mass be- 
tween the lower limit and the detection limit. Thus we in- 
tegrate the IMF between these limits and normalise. The 
IMF is generally assumed to be describable by a power law 
for supersolar masses. Here 7 is the index such that the the 
default Salpeter law gives 7 = —1.35. The parameters to be 



varied are m„ 
upper mass limit. 



the lower mass mass limit and m„ 



Pi oc 



Pi 



dm. 



the 



(7) 



(8) 



The detection limits rrii are the 84 per cent confidence 
limits. The probability of non detection if rrii exceeds rrimax 
is 1 and the probability if nii is less than mmin is 0.16. 



For the detections, the probability of the event is the 
probability that a star has this deduced mass subject to the 
errors. The distribution of the errors is unclear and the sim- 
plest way of accounting for the uncertainty is to integrate the 
IMF between the upper and lower limits set by the errors. 
However, t his skews the distrib ution towards lower masses. 
We follow ISmartt et al] (|2009l ) by instead integrating the 
IMF from the upper limit to the predicted mass and then 
in a straight line from the value of the IMF at the predicted 
mass to zero at the lower limit. If the upper error mass mi+ 
exceeds rrimax then the integral is truncated at rrimax- Simi- 
larly if the lower error mass rui- is less than rUmin then the 
integral is truncated at rrimin- 

The maximum likelihood values are of little interest 
without some measure of how probable alternatives are. To 
do this the confidence regions must be determined. For two 
parameters we have the 68, 9 and 95 pe r cent confidence 
regions when x = 2.3, 4.6, 6.2 ijPress et al.llT992ll . 

In^Cniax - ln£ = ix- (9) 

We calculated the maximum likelihood contours when 
both the upper and lower limits are varied. Initially we used 
only the detections. The non-detections favour arbitrarily 
low values for the lower limit because the slope of the IMF 
makes low-mass stars more probable. Then, with the lower 
limit fixed from the detections, the non-detections were in- 
cluded to see what effect they had on the upper limit. 

The contrast between the two models is shown in Fig- 
ure U The models without dust predicted an upper limit of 
I8+2M0 and a 95 per cent confidence limit of 25 Mq. How- 
ever, the dust models give 2lj:^M0 and, more significantly, 
a 95 per cent limit of more than 30 Mq. This means that 
we can no longer say with certainty that there is population 
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5 10 15 20 25 



'^actual 

Figure 6. Actual compared with deduced initial masses of red 
supergiants for theoretical observations in V , and K. Models 
are from the Z = 0.02 series. 



of red supergiants that do not end as Type IIP SNe. This 
upper hmit is consistent with that obtained by modelhng 
the populatio n of SN progenitor s an d accounting for binary 
evolution, i.e. lSmith et al.l(|201lh and lEldridge et all l|201lh . 

In Figure (5] we present plots of the deduced initial 
masses of the progenitors with and without dust, in order of 
increasing mass. They are contrasted with the curves repre- 
senting a population of stars following the Salpeter distribu- 
tion, with upper limits of 16.5 and 25 Mq. In both cases the 
lower hmit is 8.5 Mq. It should be noted that SN 1999ev, 
the most massive progenitor, undergoes the greatest change 
in predicted mass when dust is considered, from I8I3 
20l| /M0 . No other star is more influential in deducing the 
upper limit and our results will be more certain when we 
have more similarly massive stars. 

Recently a pro genitor to SN 2009kr was identified 
l|Fraser et al.l l2011bl ) but it is not clear yet whether it 
is a single or a binary star, a Type IIP or a Type IIL 
las- Rosa et al.l l20lrf). If it is a single star, the observa- 
tions of Mv = -7.6 ± 0.6 and Mv - Mi = 1.1 ± 0.25 imply 
initial masses of 2lli^M0 and 2315Mq from the models 
without and with dust. However, the star is more of a yel- 
low supergiant and may have evolved through the red su- 
pergiant phase (|Elias-Rosa et af]|2010l ). in which case our 
models would not be appropriate. 



4 DISCUSSION 

We have presented evidence that the red supergiant prob- 
lem is caused by aliasing of the higher masses. This is illus- 
trated in Figure [6] where we have plotted the initial mass 
of red supergiants at solar metallicity against the mass that 
would be deduced if we were to take the magnitudes of the 
dust-extincted models and compare them with the mass- 
luminosity relation from the models without dust. 

The model for circumstellar dust is fairly crude, based 
on a best fit between luminosity and dust production from 
which there is considerable deviation. It generates an aver- 
age extinction of never more than about 1 Av , whereas ob- 



serv ations have reve aled stars with several times that value 
( Massev et all |2005|[ ) . 

Detailed models of dusty winds have shown that they 
do not maintain spherical symmetry but form transient cap 
structures jWo itkc 2006). This implies that that the amount 
of extinction varies with the li ne of sight and over time. 
This variability was observed bv lMassev et al.l l|2009bl ). who 
found an average change of 0.5 mag in the V band in a 
sample of red supergiants in M31. This change occurred after 
only three years. Significantly, there was no change in the K 
band, strongly implying that the change in the magnitudes 
was driven by variable extinction. It is also likely that dust 
production varies with metallicity but we have not taken 
this into account. 

The main achievement of our dust model is to show 
that, even with an unnaturally unvarying but realistic 
amount of dust production, the red supergiant problem 
ceases to be. The increased extinction of the higher-mass 
models introduces such uncertainty into the observations as 
to make it impossible to confidently set an upper mass limit 
lower than 25-30 M©, the red supergiant upper limit. 

The upper mass limit could be determined more pre- 
cisely by obtaining more pre-explosion images in the in- 
frared, where the effect of extinction is much less. Alterna- 
tively, good spectroscopy of the detected progenitors could 
be used to calculate the total extinction. This would require 
waiting for one of the relatively small number of nearby 
well-studied red supergiants to explode. Models of the SNe 
themselves are also likely to yield more progenitor masses 
in the future and recent work in t hat direction has been 
encouraging IjPessart fc Hillieij|201lh . Until then the uncer- 
tain but significant circumstellar extinction means that we 
do not need to look for alternatives to Type IIP SNe for the 
death of red supergiants. 



ACKNOWLEDGEMENTS 

JJW and JJE would like to thank the anonymous referee for 
their suggestions. They have led to a much improved paper. 
JJW thanks the STFC for his studentship and JJE is sup- 
ported by the loA's STFC theory rolling grant. The authors 
would also like to thank Stephen Smartt and Christopher 
Tout for discussion and comments. 



REFERENCES 

Bohm-Vitense E., 1958, ZAp, 46, 108 
Cardelh J. A., Clayton G. C, Mathis J. S., 1989, ApJ, 345, 
245 

Crowther P. A., 2007, ARAA, 45, 177 

Danchi W. C, Bester M., Degiacomi C. G., Greenhill L. J., 

Townes C. H., 1994, AJ, 107, 1469 
de Jager C, Nieuwenhuijzen H., van der Hucht K. A., 1988, 

A&As, 72, 259 
Dessart L., Hillier D. J., 2011, MNRAS, 410, 1739 
Eggleton P. P., 1971, MNRAS, 151, 351 
Eldridge J. J., Izzard R. G., Tout C. A., 2008, MNRAS, 

384, 1109 

Eldridge J. J., Langer N., Tout C. A., 2011, MNRAS, 692 



© 0000 RAS, MNRAS 000, [T]-?? 




Figure 4. Maximum likeliiiood contours. Tlie daslied lines are for the detections only and the solid are when the non-detections are 
included. 



No Dust 



Dust 



1999ev - 
2003ie - 
2004dj - 
1999on - 
1999em - 
2002hh - 
2004A - 
2008bk - g 
2004am - > 
2001du - c 
20a4et - ^ 
1999g; - =1 
2007OO - 

o 

1999br -z 
2006bc - 
2006my - 
2004dg - 
2OO60V - 

2005CS 

2003gd 



15 20 
nitial mass / Mq 



1999ev 
200318 
1 999on 
1999em 
2002hh 
2004dj 
2004A 
2008bk 
2001du 
2004am 
2004et 
1999g; 
2007oa 
1999br 
2006bc 
2006my 
2004dg 
2OO60V 
2005CS 
2003gd 



15 20 
Initial moss / Mq 



Figure 5. Diagram showing the derived masses of progenitors compared to the expected mass distribution if the maximum Type IIP 
SNe progenitor mass is 16.5 M0 or 25 Mq. 



Eldridge J. J., Mattila S., Smartt S. J., 2007, MNRAS, 376, 
L52 

Eldridge J. J., Stanway E. R., 2009, MNRAS, 400, 1019 
Eldridge J. J., Tout C. A., 2004, MNRAS, 353, 87 
Elias-Rosa N., Van Dyk S. D., Li W., Miller A. A., Sil- 
verman J. M., Ganeshalingam M., Boden A. F., Kasliwal 
M. M., Vinko J., Cuillandre J.-C, Filippenko A. V., Steele 
T. N., Bloom J. S., Griffith C. V., Kleiser I. K. W., Foley 
R. J., 2010, ApJl, 714, L254 
Filippenko A. V., 1997, ARAA, 35, 309 
Eraser M., Ergon M., Eldridge J., Valenti S., Pastorello A., 
SoUerman J., Smartt S. J., Agnoletto I., Arcavi I., Benetti 



S., Botticella M.-T., Bufano F., Campillay A., Crockett 
R. M., Gal- Yam A., Kankare E., Leloudas G., Maguire K., 
Mattila S., Maund J., Salgado F., Stephens A., Tauben- 
berger S., Turatto M., 2011a, MNRAS (in press), arXiv 
1011.6558 

Eraser M., Takats K., Pastorello A., Smartt S. J., Mattila 
S., BotticeUa M.-T., Valenti S., Ergon M., SoUerman J., 
Benetti S., Bufano F., Crockett R. M., Danziger I. J., 
Maund J. R., Taubenberger S., Turatto M., 2011b, ApJ, 
714, L280 

Heger A., Fryer C. L., Woosley S. E., Langer N., Hartmann 
D. H., 2003, ApJ, 591, 288 



© 0000 RAS, MNRAS 000, [T]-?? 



10 Joseph J. Walmswell, John J. Eldridge 



Kotak R., Vink J. S., 2006, A&A, 460, L5 

Levesque E. M., Massey P., Olsen K. A. G., Plez B., Jos- 

selin E., Maeder A., Meynet G., 2005, ApJ, 628, 973 
Levesque E. M., Massey P., Olsen K. A. G., Plez B., Meynet 

G., Maeder A., 2006, ApJ, 645, 1102 
Li W., Van Dyk S. D., Filippenko A. V., Cuillandre J.-C, 

Jha S., Bloom J. S., Riess A. G., Livio M., 2006, ApJ, 

641, 1060 

Massey P., Plez B., Levesque E. M., Olsen K. A. G., Clay- 
ton G. C., Josselin E., 2005, ApJ, 634, 1286 

Massey P., Plez B., Levesque E. M., Olsen K. A. G., Silva 
D. R., Clayton G. C, 2009a, ASP, 3 

Massey P., Silva D. R., Levesque E. M., Plez B., Knut 

A. G., Geoffrey O., Clayton C, Meynet G., Maeder A., 
2009b, ApJ, 703, 420 

Mattila S., Smartt S. J., Eldridge J. J., Maund J. R., Crock- 
ett R. M., Danziger I. J., 2008, ApJl, 688, L91 

Maund J. R., Smartt S. J., Danziger L J., 2005, MNRAS, 
364, L33 

Mauron N., Josselin E., 2011, A&A, 526, A156+ 

Press W. H., Teukolsky S. A., Vetterling W. T., Flannery 

B. P., 1992, Numerical Recipes in fortran. The Art of Sci- 
entific Computing, 2nd edn. Cambridge Univ. Press 

Schrder K. P., Pols O. R., Eggleton P. P., 1997, MNRAS, 
285, 696 

Smartt S. J., 2009, ARAA, 47, 63 

Smartt S. J., Eldridge J. J., Crockett R. M., Maund J. R., 

2009, MNRAS, 395, 1409 
Smith N., Li W., Filippenko A. V., Chornock R., 2011, 

MNRAS, 412, 1522 
Stancliffe R. J., Eldridge J. J., 2009, MNRAS, 396, 1699 
Utrobin V. P., 2004, Astron. Lett., 30, 293 
Utrobin V. P., Chugai N. N., 2008, A&A, 491, 507 
— , 2009, A&A, 506, 829 

van Loon J. T., Cioni M.-R. L., Zijlstra A. A., Loup C, 

2005, A&A, 438, 273 
Vink J. S., de Koter A., Lamers H. J. G. L. M., 2001, A&A, 

369, 574 

Weidemann V., Koester D., 1982, A&A, 121, 77 
Westera P., Lejeune T., Buser R., Cuisinier F., Bruzual G., 

2002, A&A, 381, 524 
Whittet D. C. B., 2003, Dust in the Galactic Environment. 

Bristol: lOP 

WilUams K. A., Bolte M., Koester D., 2009, ApJ, 693, 355 

Woitke P., 2006, A&A, 452, 537 

Yoon S., Cantiello M., 2010, ApJl, 717, L62 



